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ABSTRACT 

Although galaxies, groups, and clusters contain ^ 10% of the baryons, many more reside in the 
photoionized and shocked-heated intergalactic medium (IGM) and in the circumgalactic medium 
(CGM). We update the baryon census in the (H I) Lya forest and warm-hot IGM (WHIM) at 10^-^ K 
traced by O VI A1032, 1038 absorption. Using robust cosmological simulations of heating, cooling, and 
metal transport, we improve the O VI baryon surveys with spatially averaged corrections for metallicity 
{Z/Zq) and O VI ionization fraction (/ovi)- Statistically, their product correlates with column 
density, {Z/ZQ)foYi w (0.015) (iVovi/lO" cm-^jO-™^ with a A^ovi-weighted mean of 0.01, which 
doubles previous estimates of WHIM baryon content. We also update the Lya forest contribution to 
baryon density out to z = 0.4, correcting for the (1 + z)^ increase in absorber density, the (1 + z)'^-'^ 
rise in photoionizing background, and cosmological proper length dlj dz. We find substantial baryon 
fractions in the photoionized Lya forest (28 ± 11%) and WHIM traced by O VI and broad-Lya 
absorbers (25 ±8%). The collapsed phase (galaxies, groups, clusters, CGM) contains 18 ±4%, leaving 
an apparent baryon shortfall of 29 ± 13%. Our simulations suggest that ~ 15% reside in hotter 
WHIM (T > 10^ K). Additional baryons could be detected in weaker Lya and O VI absorbers. 
Further progress requires higher-precision baryon surveys of weak absorbers, down to minimum column 
densities iVm > lO^^.o ^^^-2^ ^^^^ > -[^q12.5 ^^^^-2^ Noww > lO"'^ cm'^, using high-S/N data from 
high-resolution UV and X-ray spectrographs. 

Subject headings: cosmological parameters — observations — intergalactic medium — quasars: ab- 
sorption lines 



1. INTRODUCTION 

For low-redshift cosmology and galaxy formation rates, 
it is important to account for all the baryons synthesized 
in the Big Bang. Cosmologists have noted a baryon 
deficit in the low-redshift universe (Fukugita, Hogan, 
& Peebles 1998) relative to the predicted density syn- 
thesized in the Big Bang. Although this deficit could 
arise from an incomplete inventory, it could also chal- 
lenge our understanding of the thermodynamics of struc- 
ture formation and the response of the gas to accretion 
shocks and galactic outflows. Recent analysis (Komatsu 
et al. 2011) of the spectrum of acoustic peaks in the 
Cosmic Microwave Background (CMB) obtained by the 
Wilkinson Microwave Anisotropy Probe (WMAP) found 
that baryons comprise a fraction il^ = 0.0455 ± 0.0028 
of the critical matter-energy density of the universe, 
Per — (9.205 X 10^'^" g cm~'^)/iyQ, where h-jQ is the Hub- 
ble constant (i?o) in units of 70 km s~^ Mpc^^. This 
4.6% baryon fraction corresponds to a mean comoving 
density, Pf, = ribPcr — 4.24 x 10"'^^ g cm~'^, and a hydro- 
gen number density nn = 1.90 x 10~^ cm~'^, assuming a 
primordial helium mass fraction Yp = 0.2477 (Peimbert 
et al. 2007). 

Baryon inventories in the low redshift universe are un- 
certain. They are often complicated by the formation 
of galaxies and large-scale structure and by the feedback 
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from star formation in the form of ionizing radiation, 
metals, and outflows. Galaxy surveys have found ^10% 
of these baryons in collapsed objects such as galaxies, 
groups, and clusters (Salucci & Persic 1999; Bristow & 
PhiUipps 1994; Fukugita & Peebles 2004). Over the last 
15 years, substantial reservoirs of gas have been found in 
the intergalactic medium (IGM), in the halos of galax- 
ies, and in the circumgalactic medium (CGM). There 
are often semantic problems in defining the CGM, as gas 
blown out of galaxies or material within the virial ra- 
dius (Tumlinson et al. 2011; Prochaska et al. 2011). Of 
the remaining 80-90% of cosmological baryons, approxi- 
mately half can be accounted for in the low-z IGM (Shull 
2003; Bregman 2007; Danforth & Shull 2008; Danforth 
2009) including the warm-hot IGM (or WHIM). Uhravi- 
olet spectroscopic surveys of Lya and O VI have identi- 
fied substantial numbers of absorbers (Danforth & Shull 
2008; Tripp et al. 2008; Thom & Chen 2008), but claimed 
detections of hotter in X-ray absorption by O VII (Nicas- 
tro et al. 2005) remain controversial (Kaastra et al. 2006; 
Yao et al. 2012). Unfortunately, X-ray spectra still have 
not confirmed the potential large reservoir of baryons at 
T > 10^ K, suggested by cosmological simulations. 

An inefficient distribution of collapsed baryons vs. dis- 
tributed matter is a prediction of nearly all cosmolog- 
ical simulations of large-scale structure formation (Cen 
& Ostriker 1999, 2006; Dave et al. 1999, 2001; Smith 
et al. 2011; Tepper-Garcia et al. 2011). These N-body 
hydrodynamical simulations suggest that 10-20% of the 
baryons reside in collapsed objects and dense filaments, 
with the remaining 80% distributed over a wide range of 
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phases in baryon overdensity (Ab — Pbl~Ph) ^^^d temper- 
ature (T). Thermodynamic considerations suggest that 
shock-heated WHIM at 2; < 1 is a natural consequence 
of gravitational instability in a dark-matter dominated 
universe. This hot gas is augmented by galactic-wind 
shocks and virialization in galaxy halos. 

In this paper, we improve the analysis of baryon 
content in both photoionized and shocked-heated IGM 
phases and assess the accuracy of observational and theo- 
retical estimates. Because Lya absorption surveys probe 
the neutral component of diffuse, photoionized filaments, 
assessing the baryon content requires a correction for the 
neutral fraction, /hi- We derive more accurate photoion- 
ization corrections based on the (1 -f z)*-^ rise in the 
metagalactic ionizing background out to 2: « 0.7, as well 
as the (1 -f z)^ increase in absorber density and changes 
in cosmological proper length, dljdz. We also use cos- 
mological simulations to correct the O VI surveys for ion- 
ization fraction (/ovi) and metallicity {Z/Zq). Finally, 
we use our group's critically evaluated catalog of H I and 
O VI absorption lines (Tilton et al. 2012) from our HST 
Legacy Archive project on IGM data taken by the Hubble 
Space Telescope (HST) and Far Ultravolet Spectroscopic 
Explorer (FUSE). 

To constrain the distribution of Z and /ovii we 
use cosmological hydrodynamic simulations of IGM 
heating, cooling, and metal transport to find the 
column-density weighted average for their product, 
(Z/Zq) /ovi- Whereas previous work assumed con- 
stant values, Z/Zq = 0.1 and /ovi = 0-2, our new 
simulations provide spatially averaged corrections for 
metallicity, O VI ionization fraction, and covariance of 
IGM parameters (T, p, Z). Our computed mean value 
{{Z/Zq) f ovi) = 0.01 is half the previously assumed 
value, thereby doubling previous estimates of the baryon 
census traced by O VI. In Section 2 we describe the 
simulations and their results for the O VI distribution 
in column density, gas temperature, baryon overdensity, 
metallicity, and ionization fraction. In Section 3, we as- 
sess the corrections for ionization and metallicity, ap- 
plied to the O VI and Lya absorbers, and we derive val- 
ues of from recent surveys of IGM and CGM phases. 
Variations in these factors are produced by WHIM ther- 
modynamics, gas temperature (T), baryon overdensity 
(Ab), and O VI column density (iVovi)- In Section 4 
and two Appendices, we summarize the current baryon 
census with uncertainties on each component. 

2. COSMOLOGICAL SIMULATIONS OF THE WHIM 

The primary simulation analyzed in this work is run 
50_1024_2 of Smith et al. (2011) with "distributed feed- 
back" and performed with the adaptive-mesh-refinement 
-|- N-body code Enzo (Bryan & Norman 1997; O'Shea 
et al. 2004). To check convergence and robustness of 
our results, we also looked at other runs, particularly 
run 50_1024_1 with "local feedback" . Post-processing of 
the simulations was carried out using the data analysis 
and visualization package ytj^ documented by Turk et al. 
(2011). The simulation has a box size of 50h^^ Mpc co- 
moving, with 1024'^ grid cells and dark matter particles, 
giving it a dark-matter mass resolution of 7 x 10^ Mq and 
spatial cells of 50h~^ kpc. Radiative cooling is included 
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by solving for the non-equilibrium chemistry and cooling 
of atomic H and He, coupled to tabulated metal cool- 
ing rates computed as a function of density, metallicity, 
temperature, electron fraction, and redshift in the pres- 
ence of an ionizing metagalactic radiation background. 
To mimic the effects of reionization, we included a spa- 
tially uniform, but redshift-dependent radiation back- 
ground given by Haardt & Madau (2001). The infiuence 
of the radiation background through photoheating and 
photoionization is included in the cooling of both the 
primordial and metal species. 

Our Enzo unigrid simulations (Smith et al. 2011), post- 
processed in the current paper, are among the best cur- 
rent work in describing the temperature, metallicity, and 
ionization state of the hot gas (WHIM) and photoion- 
ized gas. Our results have been checked for convergence 
(Smith et al. 2011), and they are robust. We have tested 
the results from several runs with different box sizes, res- 
olutions, and modes of feedback (local and distributed), 
finding similar results, as described later. Our use of un- 
igrid (rather than adaptive-mesh-refinement) hydrody- 
namics was deliberate, to provide optimal hydrodynamic 
resolution of the IGM and WHIM. It is well estabhshed 
that SPH (smooth particle hydrodynamic) codes provide 
poor hydrodynamic resolution of shocks and instabili- 
ties. Agertz et al. (2007) noted that Eulerian grid-based 
methods are able to resolve and treat important dynam- 
ical instabilities such as Kelvin-Helmholtz or Rayleigh- 
Taylor, whereas these processes are poorly resolved by 
SPH techniques. Vazza et al. (2007) pointed to the sig- 
nificantly different phase diagrams of shocked cells in grid 
codes compared to SPH, with sizable differences in the 
morphologies of accretion shocks between grid and SPH 
methods. 

We have employed a sophisticated treatment of metal 
cooling and heating (Smith, Sigurdsson, & Abel 2008; 
Smith et al. 2011). The SPH models by Dave & Oppen- 
heimer (2009) and Dave et al. (2011) find IGM phases 
with markedly different temperatures and metallicities, 
but we believe these results arise from their metal- 
injection schemes and the lack of mixing. None of the 
other major simulation groups (Gen & Ostriker 2006; 
Gen & Chisari 2011; Tepper-Garcia et al. 2011; Smith 
et al. 2011; Gen 2012) can reproduce their results on 
IGM temperature and ionization state. Our 1024'^ simu- 
lations have better hydrodynamic resolution than all of 
the Tepper-Garcia et al. simulations (512"^) and Dave & 
Oppenheimer models (384'^). Our paper was the only one 
to test for convergence, and we concluded that our 1024^^ 
models were robust. It has sometimes been suggested 
that the use of "unigrid" Enzo calculations casts doubt 
on the validity of the results for metallicity production 
and transport, compared to adaptive- mesh-refinement 
(AMR). Just the opposite is true. In fact, AMR does 
not provide a better computation of the global star for- 
mation rate, as we examined (Smith et al. 2011). In 
these simulations, star formation seems to be dependent 
on the force resolution on the root grid. Stars do not 
form in the lower-resolution simulations because the ha- 
los do not collapse correctly at early times when there is 
no refinement occurring in the box. By the time refine- 
ment happens, it is already too late; the collapse of the 
halos has already been delayed, owing to the poor force 
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resolution. Wc have verified this by running simulations, 
with and without AMR, at the same resolution on the 
top grid. We find almost the exact same star-formation 
history. For this reason, we decided to use high reso- 
lution everywhere and put our computational resources 
into a large unigrid: 1024^ in this paper, and moving to 
1536-'' in recent work at z « 6 - 10 (ShuU ct al. 2012a). 

We use a modified version of the star-formation rou- 
tine of Cen & Ostriker (1992). Star particles, represent- 
ing the combined presence of a few million solar masses 
of stars, are formed when the following three conditions 
are met: the total density of a grid cell is above a cer- 
tain threshold, a convergent flow exists (negative velocity 
divergence) , and the cooling time is less than the dynam- 
ical time. Star particles return feedback to the grid in 
the form of metal-enriched gas and thermal energy. We 
use a distributed-feedback method (Smith et al. 2011) 
in which material and energy are distributed over a 27- 
cell cube (three cells on a side) centered on the location 
of the star particle. This is done to mitigate the over- 
cooling that occurs when feedback is deposited into a 
single grid cell, which raises the temperature, gas den- 
sity, and cooling rate to unphysically high values. Over- 
cooling also causes the winds that should be transporting 
metals into the IGM to fizzle out and remain confined to 
their sites of origin. Smith et al. (2011) showed that this 
model is able simultaneously to provide good matches to 
the global star formation history and the observed num- 
ber density per unit redshift of O VI absorbers (Dan- 
forth & ShuU 2008). If all three above conditions are 
satisfied, a "star particle," representing a large collec- 
tion of stars, is formed within the grid cell with a total 
mass, TO* = /*mceii(Ai/tdyn)- Here, /* « 0.1 is the 
star-formation efficiency, TOceii is the baryon mass in the 
cell, tdyn is the dynamical time, and At is the hydrody- 
namical timcstcp. Subsequently, this much mass is also 
removed from the grid cell as the star particle is formed, 
ensuring mass conservation. Although the star particle is 
formed instantaneously within the simulation, feedback 
is assumed to occur over a longer time scale, which more 
accurately reflects the gradual process of star formation. 
Stellar feedback is represented by the injection of ther- 
mal energy and the return of gas and metals to the grid, 
in amounts proportional to Awsf. We assume that a 
fraction (25%) of the stellar mass is returned to the grid 
as gas. The thermal energy and metals returned to the 
grid are e = e(ATOsf)c^ and TOmctais = {Arnsi)y, where 
e « 10~^ is the ratio of rest-mass energy to thermal en- 
ergy and y « 0.025 is the metal yield. The rate of star 
formation and rate at which thermal energy and metals 
are injected into the grid by the star particle peak after 
one dynamical time, then decays exponentially. 

In the current project, wc seek to understand the ther- 
mal and ionization state of the IGM, including the dis- 
tribution and covariance of metallicity, temperature, and 
O VI ionization fraction. In our simulations, the physi- 
cal properties of O VI absorbers and the degree to which 
they trace the WHIM are in good agreement with other 
recent simulations (Tepper-Garcia et al. 2011; Cen & 
Chisari 2011). However, they differ significantly from 
those of Oppenheimcr & Dave (2009, 2011), who found 
that O VI and Ne VIII originate almost exclusively in 
warm (T « 10^ K) photoionized gas. Tepper-Garcia 



ct al. (2011) suggested that much of this difference arises 
because Oppenheimcr & Dave (2009) neglected the effect 
of photoionization on metal-line cooling. Oppenheimcr 
et al. (2012) investigated this claim by running additional 
simulations in which photoionization reduced the metal 
cooling. This reduction resulted in a small, but insuffi- 
cient number of O VI absorbers associated with WHIM 
gas. Instead, they point to the fact that they do not in- 
clude the mixing of metals from their galactic- wind par- 
ticles, allowing feedback to take the form of cold, heavily 
enriched clouds. 

In examining the various interpretations of the high 
ions (O VI and Nc VIII), we have identified several key 
differences between the two codes: the Oppenhcimer- 
Dave smooth particle hydrodynamic (SPH) code and our 
grid-based approach with Enzo. The disagreement be- 
tween the IGM temperatures and ionization mechanisms 
appear to arise from four effects: (1) Different meth- 
ods of injecting energy and metals (mixed or unmixed); 
(2) Over-cooling of unmixed, metal-enriched gas at high 
density and high metallicity; (3) Differences in shock cap- 
turing (and shock-heating) between SPH and grid codes; 
(4) Differences in photoionization rates through the as- 
sumed radiation fields at hiy = 100 — 250 eV. All four 
possibilities merit careful comparative studies, although 
we believe the resolution of this O VI controversy may 
ultimately hinge on understanding the nature of galactic 
winds and their ability to mix heavy elements into IGM. 

3. CENSUS OF BARYONS IN DIFFERENT THERMAL 
PHASES 

The "Lyman-a forest" of absorption lines appears to 
contain ~30% of the low- 2: baryons (Penton et al. 2000, 
2004: Lehncr et al. 2007). Another 30 40% is predicted 
by simulations to reside in shock-heated gas at 10^ K to 
10^ K (WHIM). Owing to its low density, the WHIM 
is difficult to detect in emission (Soltan 2006). More 
promising are absorption-line studies, using the high ion- 
ization states of abundant heavy elements with resonance 
lines in the far- ultraviolet (C IV, N V, O VI), extreme ul- 
traviolet (O IV, O V, Ne VIII), and soft X-ray (O VII, 
O VIII, N VI, Ne IX). The low-z WHIM has most effec- 
tively been surveyed in the O VI lines at 1031.926 A and 
1037.617 A (Daiiforth & ShuU 2005, 2008; Tripp ct al. 
2008; Thom & Chen 2008), which probe the temperature 
range l0^-^~^-^ K in coUisionally ionized gas. Danforth 
& ShuU (2008) measured the column densities of 83 O VI 
absorbers and estimated that 8.6 ± 0.8% of the baryons 
reside in this phase, assuming constant correction fac- 
tors for the metallicity (Z « 0.1 Zq) and O VI ionization 
fraction (/ovi = 0.2). A few detections of Ne VIII have 
also been reported (Savage et al. 2005; Narayanan et al. 
2009, 2011; Meiring et al. 2012) probing somewhat hotter 
gas. Weak X-ray absorption lines are difficult to detect 
with the current throughput and spectral resolution of 
spectrographs on Chandra and XMM/Newton. Possible 
X-ray detections of hotter gas at (1 — 3) x 10^ K have 
been claimed, using absorption lines of helium-like O VII 
A21.602 (Nicastro et al. 2005a,b, 2008; Buote ct al. 2009; 
Fang et al. 2010; Zappacosta et al. 2010) and hydro- 
genic O VIII A18.969 (Fang et al. 2002, 2007). Most of 
these Chandra detections remain controversial and un- 
confirmed by the XMM-Newton satellite (Kaastra et al. 
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Fig. 1. — Distribution of IGM in temperature T and baryon overdensity A;, = pi,/pfj, color-coded by baryon mass fraction. This 
distribution shows the same thermal phases seen in simulations by other groups and commonly labeled as warm (diffuse photoionized gas), 
WHIM, and condensed. 



2006; Williams et al. 2006; Rasmussen et al. 2007). Re- 
cent analyses of spectroscopic data on Mrk 421 fails to 
detect any WHIM gas at the claimed redshifts {z — 0.01 
and 0.027), either in broad Lya absorption (Danforth 
et al. 2011) from high-S/N data from the Cosmic Ori- 
gins Spectrograph (COS) on the Hubble Space Telescope 
(hST) or in O VII (Yao et al. 2012) in Chandra data. 

The post-processed results of our simulations relevant 
to O VI are illustrated in Figures 1-7. Figure 1 shows 
the temperature-density phase diagram of the multiphase 
IGM, color-coded by baryon mass fraction. The com- 
monly found features in (T, A^) plots include the dif- 
fuse Lya absorbers (T = lO^ o-^-^ k and = 10"^ to 
10+^-^), a condensed phase (Ab > lO'^), and a shocked- 
heated plume of WHIM (T = 10^"^ K and A,, = 10° to 
10^-^). The gas traced by O VI resides at 5 < logT < 6, 
marked in Figure 2 for the cumulative distribution of 
mass and metals vs. temperature in our simulations. The 
two curves show results from simulations with distributed 
feedback (energy injected into 27 cells) and local feed- 



back (single cell). The similarity of the mass distribu- 
tion curves is evidence of model robustness. Figure 3 
shows the cumulative mass distribution of O VI versus 
column density, indicating that a significant fraction of 
O VI is traced by weak absorbers with logiVovi < 13-5 
(equivalent width Wx < 40 mA for stronger 1032 A line). 

Our simulations show a large range of O VI proper- 
ties, including metallicity, temperature, and ionization 
fraction. Figures 4 and 5 show the range and covariance 
of the two correction factors: metallicity in solar units 
(Z/Zq) and O VI ionization fraction /ovi- The ioniza- 
tion fraction includes the effects of both collisional ion- 
ization and photoionization, as discussed by Smith et al. 
(2011). Figure 4 color-codes these distributions in A^ovi, 
while Figure 5 illustrates their distribution in tempera- 
ture and baryon overdensity. The wide range of physical 
conditions in which O VI exists demonstrates the multi- 
phase character of this ionized gas. The mechanisms that 
produce O VI (113.87 eV is needed to ionize O V) in- 
clude collisional ionization at log T > 5 and photoioniza- 
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tion by the hard (EUV) metagalactic background. The 
dashed hne in Figure 5 shows the locus of points (log T 
and logAb) at which coUisional ionization equals pho- 
toionization. This calculation is based on O V collisional 
ionization rates from ShuU & Van Steenberg (1982). The 
O V photoionization rate was derived, assuming a cross 
section cr(i;) = (0.36x10-^^ cin'^){E/lU.87 eY)'^-^ and 
the EUV radiation field at 8-10 ryd from Figure 13 of 
Haardt & Madau (2012). 

Figure 6 shows the distribution of the product, 
(Z/Zq) /ovI: throughout the simulation. Color-coded by 
O VI mass fraction, this plot exhibits bimodality of this 
product with temperature and baryon overdensity. In the 
deepest-red portions, representing high O VI mass frac- 
tion, one finds regions with large values, (Z/Zq) /ovi ~ 
lO""'^'^, in high-density filaments with A{, ^ 100, and 
other regions with lower values, {Z/Zq) /ovi ~ 10~^'^ at 
Afe f« 10. In converting the column densities of O VI ab- 
sorbers to baryons, one therefore should not adopt single 
corrections for metallicity and ionization fraction. How- 
ever, as we discuss in Section 3.1, one can apply statis- 
tical corrections that correlate the product {Z/Zq) /ovi 
with O VI column density (Figure 7). Integrating over 
the O VI column-density distribution, we find that this 
product is a factor of two lower than the previously as- 
sumed value of 0.02. This doubles the fraction of WHIM 
baryons traced by O VI from 8-9% to 17 ± 4%. 



3.1. 



Warm-Hot IGM Probed in VI Absorption 



As noted, the metallicity and O VI ionization fraction 
vary throughout the grid in the simulated IGM. Figures 6 
and 7 show that the product {Z/Zq) f ovi correlates with 
Novii with weaker absorbers having systematically lower 
values. This is primarily an effect of spatial variations in 
the metallicity, but also variations in the O VI ionization 
fraction produced by spatial fluctuations in WHIM tem- 
perature and contributions from photoionizing radiation 
on low-density IGM. The baryon content in O Vl-traced 
WHIM is given by an integral over O VI column den- 
sity (iV), now with the corrections for metallicity and 
ionization fraction placed inside the integral: 



(OVI) 



cp„(0/H), 



(W{N) 
dz 



N dN 



a. Zo{N)foYl{N) 

Here, /i;, = 1.33m// is the mean baryon mass per hydro- 
gen nucleus, accounting for helium. Earlier surveys that 
used O VI as a baryon tracer assumed constant values of 
the ionization fraction, /ovi ~ 0-2 (its maximum value 
in collisional ionization equilibrium at logTmax = 5.45) 
and metallicity, Z/Zq =0.1, relative to the solar oxygen 
abundance, (O/H)© = 4.90 x 10"'' (Asplund et al. 2009). 

From our simulations (Figure 7) , the product of metal- 
licity and O VI ionization fraction has a statistical power- 
law dependence on O VI column density, scaling as iV", 



{Z/Zq) fowl = (0.015) (iVovi/10' 
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Integrating over the distribution dM /dz oc N ^ , we find 
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Fig. 2. — Cumulative distributions of IGM mass (top panel) and 
metals (bottom panel) vs. temperature T from our simulations, run 
with both "local and distributed" feedback (Smith et al. 2011). 
We only show IGM baryons with overdensities A(, < 1000; the 
remaining mass and metals are in galaxies and collapsed phase. 
Vertical dashed lines mark the temperature range (10^ K to 10^ K) 
probed by O VI and Ne VIII, and horizontal lines show intercepts 
for the two distributions. 

For 7 « 0.70 and the observed /3 sa 2.0 (Danforth & 
ShuU 2008), the integral increases at the low end as 
^min''- "^^^ column-density weighted mean of the prod- 
uct /ovi {Z/Zq) = 0.01, a factor of two smaller than 
previously assumed. For our standard integration range, 
13.0 < logiVovi ^ 15.0, this correction doubles the num- 
ber of baryons in the O Vl-traced WHIM, compared to 
previous assumptions. As noted earlier, we found similar 
results in all our simulations. 

Previous O VI surveys and baryon estimates were 
quantified by an absorption-line frequency, cW /c?z, per 
unit redshift and an O Vl-traced baryon fraction, 

r2[°^^Vf^b- Using 40 O VI absorbers seen with FUSE, 
Danforth & ShuU (2005) found dJ\f/dz w 17 ± 3 for col- 
umn densities 13.0 < Nowi < 14.5. The 10^^ cm'^ 
lower limit corresponds to 12.5 mA equivalent width 
in O VI A1032. The 2005 census gave an O Vl-traced 
baryon fraction of at least 4.8 ± 0.9% (statistical error 
only) for a correction factor /qvi {Z/Zq) = 0.02. Dan- 
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Fig. 3. — Cumulative distribution of O VI by mass from our 
WHIM simulation, for column densities A^ovi = 10^^ to 10^^ 
cm"'^. Solid black line shows cumulative mass fraction, labeled 
on left vertical axis. Gray histogram shows cumulative number of 
O VI absorbers, labeled on right vertical axis. Red points show ob- 
served cumulative column-density distribution (Danforth & ShuU 
2008) labeled on right vertical scale. A significant fraction of O VI- 
traced WHIM resides in weak absorbers, although the distribution 
flattens at logAfovi < 13.5. 
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Fig. 4.— Distribution of IGM metallicity (Z/Zq) and O VI ion- 
ization fraction (/ovi) color-coded by O VI column density. Note 
the wide range and covariance of individual factors, whose prod- 
uct, /ovi (Z/Zq), correlates with O VI column density (color bar 
along right). 



forth & ShuU (2008) used HST/STIS data on 83 O VI 
absorbers to find dAf/dz « 15.0^2 q integrated to 30 
mA, with dAf /dz sa 40tg'' integrated to 10 niA. Their 
derived WHIM baryon fractions were 7.3 ± 0.8% and 
8.6 ± 0.8%, integrated to 30 mA and 10 mA, respec- 
tively. Tripp et al. (2008) found a similar fine frequency, 
dM jdz K, 15.612 4 51 intervening O VI absorption 
systems integrated to 30 mA, while Thom & Chen (2008) 
found dM jdz sa 10.4 ± 2.2 for 27 O VI absorbers down to 



30 mA. The latter survey was more conservative, requir- 
ing detection of both the stronger O VI line at 1032 A 
and the weaker line at 1038 A. They also had fewer sight 
lines (16) and a much smaller number of O VI absorbers 
(27) in their survey. Their line frequency is smaller than 
those in other O VI surveys, perhaps because of small- 
number statistics or owing to the two-line requirement. 

As part of a Hubble Archive Legacy project (PI: Shull, 
HST-AR-11773.01-A), the Colorado group has reana- 
lyzed T/ST/STIS data on IGM absorption fines (Tilton 
et al. 2012), finding 746 H I absorbers and 111 O VI ab- 
sorbers. Through a critical evaluation of data in the lit- 
erature and comparison to high-S/N COS spectra when 
available, we corrected line identification errors in Dan- 
forth & ShuU (2008) and other surveys. The O VI line 
frequency from TiUon et al. (2012) is dM jdz w 22.2^2:4 

integrated to 30 mA, with a slope /3 = 2.075±0.119. This 
new survey provides more reliable baryon fractions for 
the O Vl-traced WHIM of 7.2±0.8% (integrated down to 
30 mA equivalent width) and 8.6 ±0.7% (to 10 mA), us- 
ing the old correction factor, /ovi {Z/Zq) — 0.02. Taken 
as a whole, these O VI surveys suggest a WHIM baryon 
fraction of 8-9%, assuming the old values for metallicity 
and ionization fraction. With our new correction fac- 
tors, /ovi {Z/Zq) = 0.01, the baryon fractions double. 
We therefore adopt an O Vl-traced baryon fraction of 
17 ± 4%, where we have increased the error to account 

for systematic uncertainties. The value of is dom- 

inated by weak absorbers with TVovi < lO^'^'^ cm~^, 
where the absorption-line statistics become uncertain. 
An accurate O VI census will require measuring even 
weaker O VI A1032 absorption lines, with column densi- 
ties below 10^'^ cm~^, corresponding to equivalent widths 
Wx = (12.5 mA)(A^ovi/10^^ cm~^). Deep spectroscopic 
surveys in O VI can also ascertain where the distribution 
of O VI absorbers flattens (Figure 3). 

3.2. Diffuse Photoionized Lya Filaments 

In this subsection, we provide details on the assump- 
tions, described briefly in Penton et al. (2000), for the 
photoionization corrections to the H I (Lya) absorbers. 
In photoionization equilibrium, the density of neutral hy- 
drogen in low-density gas depends on the ionizing back- 
ground, gas density, and electron temperature, all of 
which evolve with redshift and have spatial fluctuations 
throughout the IGM: 



"-HI 



He riH ot-u'^iT) 



H 



(4) 



Here, rig = (1 -I- 2y) = 1.165n// is the electron density 
for fully ionized gas, nu is the total density of hydrogen 
nuclei, and y = nHc/nn = [(>"/4)/(l - Y)] « 0.0823 
is the helium-to-hydrogen ratio by number, assuming 
Y — 0.2477 helium abundance by mass. For the low 
values of uh and Nhi in IGM absorbers, we adopt the 
hydrogen case- A radiative recombination rate coefficient. 



{A) 



(T) w (2.51 X 10- 



cm^ s-i)r. 



4 3 '^^^, scaled to an 
4.3 -w\rp^ ^ characteristic 



electron temperature T = (10 K)T4 3 
of low-metaUicity IGM (Donahue & ShuU 1991). The 
H I photoionization rate depends on the metagalactic 
radiation field, with specific intensity Ii, = /o(j^/i'o)~"° 
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Fig. 5. — Distribution of IGM temperature versus baryon overdensity, A;, = p^/pj, color-coded by O VI mass fraction. In this phase 
space, we identify the WHIM (T > 10^ K), warm, diffuse photoionized gas (T < 10^ K and Aj < 1000), and condensed gas (A;, > 1000). 
Dashed line shows locus at which coUisional ionization equals photoionization (coUisional ionization dominates above the dashed curve). 



referenced to the Lyman limit, hi^o = 13.6 eV. Here, 
(as) « 1.6 — 1.8 is the mean QSO spectral index be- 
tween 1.0-1.5 ryd (Teller et al. 2002; Shuh et al. 2012b). 
The frequency-integrated photoionization rate is given 
by the approximate formula, ~ [47r/o(To/^(cts + 3)] ~ 
(2.49 X 10-" s-i)/_23(4.8/a^ + 3), where Iq = (lO'^^ 
erg cm^^ s^^ Hz^^ sr^^)/_23- Hereafter, we com- 
bine the two parameters, Iq and as, into a single scal- 
ing parameter, r_i4 = (F/z/lO-^'' s"^), for the hy- 
drogen photoionization rate. Previous calculations of 
the low-2 ionizing intensity (ShuU et al. 1999) found 
w 3.2+^;^ X 10"" s^i and noted that low-z obser- 
vational constraints were consistent with values in this 
range. More recent calculations (Haardt & Madau 2012) 
are consistent with this value and other estimates at 
z « 0. Figure 8 illustrates current estimates of the inten- 
sity of the metagalactic ionizing radiation field starting 
at 1 ryd (H I photoionization edge) , continuing to 4 ryd 



(He II edge) and beyond, including ionization potentials 
needed to produce some of the higher metal ions, such 
as C IV (47.87 eV = 3.52 ryd), O VI (113.87 eV = 8.37 
ryd), O VII (138.08 eV = 10.15 ryd), Ne VIII (207.2 eV 
= 15.24 ryd), and O VIII (739.11 eV = 54.35 ryd). 

The Lya absorbers are thought to arise as fluctuations 
in dark-matter confined clumps or filaments, which we 
approximate as singular isothermal spheres with density 
profiles, n}i{r) — no{r /ro)~'^ normalized to a fiducial 
radius tq. For a sight line passing through an absorber at 
impact parameter p, the total baryon mass within radius 
r = p is Mb{p) = 47r(Ubno''oP: where /ifc = (1 -I- Ay)mH ~ 
LSSniH is the mean baryon mass per hydrogen. The 
H I column density can be derived, using the density- 
squared dependence of nni and integrating through the 
cloud along pathlength £ at impact parameter p, where 
£2 _ j,2 _ substitute p = r COS (p and £ = ptant/) 
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Fig. 6. — Distribution of the metallicity-ionization product, 
{Z / Zq) fo\i, versus baryon overdensity, A;, = pij/ji^, color-coded 
by O VI mass fraction. The broad distribution of the prod- 
uct, from 0.001 to 0.1 has local enhancements (deep red) in low- 
metallicity regions (overdensities A;, Ki 10) and high-metaUicity 
regions (A;, 100). The column density weighted mean of this 
product is 0.01 (Figure 7). 

for the angle (j) between the directions of p and r, 

-4 



2 / n, 





dl 



(5) 



Solving for the quantity tioTq, we can express the ab- 
sorber mass within r = p as 



Af6(p) = 47r^hP^/2 



2rgiVHi(p) 
7r(l + 2y)a^^^ 



1/2 



:(1.09 X l(f Mq) 



l/2pl/2 5/2^0.363 



14 



-14 



^'l00-'4.3 



(6) 



In the above formula, we have scaled the impact param- 
eter p = (100 kpc)pioo using a 100-kpc characteristic 
scale length of Lya absorbers at fiducial column density 
Nm = (10" cm-2)iVi4; see Penton et al. (2000) and 
references therein. 

Next, we calculate the baryon density in the Lya ab- 
sorbers as a fraction of the cosmological closure density, 
Per = (SHq/SttG) at z = 0. Because our Lya survey 
extends to z = 0.4, we include corrections for cosmo- 
logical evolution in the space density of absorbers, (j>{z), 
and the hydrogen photoionization rate, Tniz). We begin 
with the standard expression for the number of absorbers 
per unit redshift. 



dz 



{l + z)H{z) 



7r[p(z)]2 0(z) 



(7) 



where H{z) = Ha[^„,{l + zf + f^A]^/^ is the Hubble 
parameter at redshift z in a flat ACDM cosmology, and 



the absorber space density is ^(z) = 4>o{\ -\- z)^ . The 
absorption-line frequency and impact parameter p cor- 
respond to the H I column density and absorber mass 
given in Equations (5) and (6). The baryon mass den- 
sity in Lya absorbers at redshift z is written as the prod- 
uct of absorber mass, Af(,, times the absorber space den- 
sity, 4>{z)^ integrated over the distribution in H I col- 

umn density. We define the closure parameter, fij^ — 
Pb(0)/Pcr(0) at 2 = 0, where pb{z) = ph(0)(l + zf. As- 
sembling all the terms, we find a closure parameter in 
Lya absorbers of column density A^hi, 



P,(Hi) _ (t^gMb 



Per 

32(2^)1/2 



cH{z){l + zf 



Th{z)p{z)N, 



-,1/2 



HI 



(1 



(8) 



Recent calculations of the metagalactic ionizing back- 
ground (Haardt & Madau 2012) show that the hydrogen 
photoionization rate rises rapidly from redshifts z = Q 
to z = 0.7. We have fitted the rate to the convenient 
formula Th{z) = (2.28 x IQ-^^ s-i)(l + zf-'^. We have 
little data on the redshift evolution of the characteristic 
absorber scale length, p, other than theoretical expecta- 
tions for the gravitational instability of filaments in the 
cosmic web. In the following calculation, we assume that 
p = (100 kpc)pioo remains constant with redshift. We 
rewrite Equation (8) as an integral over column density. 



(HI). 



(9.0 X 10" 



[r!„(l + z)3 + r!A]V2 



N„ 



N„ 



dNilogNni) 
dz 



Nli^ d{\og Nm) . (9) 



The ionization rate, L/f (z) cx (1 + z)**'^, enters Equation 
(8) as the square-root, almost exactly compensating for 
the (1 -I- z)2 factor in the denominator. 

From UV spectrographic surveys of intergalactic Lya 
absorbers, we now have a reasonable understanding of 
the distribution of H I column densities in the diffuse 
Lya forest (Penton et al. 2000, 2004; Danforth & Shull 
2008; Lehner et al. 2007; Tilton et al. 2012). Penton 
et al. (2004) applied photoionization corrections to a 
survey of 187 Lya absorbers over redshift pathlength 
Az — 1.157 and found a distribution in H I column den- 
sity f{N) cx N-^ with P = 1.65 ± 0.07 and a 29 ± 4% 
contribution to the total baryon content. Danforth & 
Shull (2008) surveyed 650 Lya absorbers with HST/STIS 
over the range 12.5 < logA'ni < 16.5 with total path- 
length Az = 5.27. Their H I distribution was similar, 
with (3 — 1.73 ± 0.04 and a fractional contribution of 
28.7 ± 3.7% (statistical errors only) to the baryon cen- 
sus. Although the formulae for baryon fractions were 
stated correctly (Equations 7-10 in Danforth & Shull 
2008), the values listed in their Tables 12 and 13 were 
computed with incorrect cosmological corrections and 
overestimated values of fif, for H I and O VI. Our new 
T/ST archive survey (Tilton et al. 2012) properly includes 
these effects, as well as the redshift evolution of Th{z), 
to find /3 = 1.68 ± 0.03 for 746 Lya systems with column 
densities between 12.5 < logA^ni < 16.5 over pathlength 
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Fig. 7. — Values of the product, (Z/Zq) foyi, for the WHIM simulation with distributed feedback overlaid with a power-law fit (solid 
red line) to the ensemble, (0.015)[A^ovi/10'^^ cm"-^]'^-^'^. Results from simulations with local feedback are similar. Integrating over the 
observed distribution in A^ovii find a weighted value of O.Of. Horizontal (blue) line shows previously assumed value of 0.02 for this 
product (Z/Zq = 0.1 and /ovi = 0.2). Typical variance (Icr) about this fit is 0.3-0.4 dex, shown by red dashed lines. 
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Az ~ 5.38. Their new derivation of fi^ is consistent 
with 24-30% of the baryons residing in the Lya forest and 
partial Lyman-Umit systems. For the current census, we 

f HI) 

adopt a baryon fraction of fij^ ^ = 28 ± 11%. The latter 
error bars include systematic effects, discussed in greater 
detail in Appendix A. 

4. BARYON CENSUS AND FUTURE SURVEYS 

We now summarize the current status of the low-z 
baryon census and discuss requirements for future sur- 
veys. Details of the individual baryon contributions are 
described in Appendix B, giving previous estimates and 
uncertainties. Figure 9 shows a pie chart of the current 
observable distribution of low-redshift baryons in vari- 
ous forms, from collapsed structures to various phases of 
the IGM, CGM, and WHIM. These slices show the con- 
tributions, il^*' /r2[*°*\ to the total baryon content from 
components (i). Measurements of Lya, O VI, and broad 
Lya absorbers, together with more careful corrections for 
metallicity and ionization fraction, can now account for 
~ 60% of the baryons in the IGM. An additional 5% may 
reside in circumgalactic gas, 7% in galaxies, and 4% clus- 
ters. This still leaves a substantial fraction, 29 ± 13%, 



unaccounted for. We have assigned realistic errors on 
each of the "slices of the baryon pie" , most of which in- 
volve systematic uncertainties in the parameters needed 
for the ionization corrections, metallicity, and geometric 
factors (cloud size). It is possible that the baryon inven- 
tory could change as a result of better determinations of 
these parameters. However, most numerical simulations 
including ours (Figure 2) suggest that a substantial reser- 
voir (~ 15%) of hot baryons exists in the hotter WHIM 
(T > 10*^ K). 

What observations and theoretical work are needed 
to make progress on error bars or baryon detections? 
First, we need more precise UV absorption-line surveys 
to measure O VI and Lya absorbers to lower column 
densities. As described in Section 3, the numbers of ab- 
sorbers in current surveys become increasingly uncertain 
at column densities log A^hi < 13.0 and log Aqvi < 13.5. 
The current surveys integrate below these levels, but 
our experience using COS to re-examine earlier Lya and 
O VI detections with FUSE, GHRS, and STIS, suggests 
that some of these weak absorbers are unconfirmed with 
high-S/N data. In our HST Archive Legacy Project, 
Tilton et al. (2012) reanalyzed O VI and Lya data from 
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Fig. 8. — Compilation of the mean AGN spectral energy distri- 
bution (SED), pinned at 1 ryd, with /q = 2 X 10~^^ erg cm~^ s"^ 
Hz~^ sr~^ and two different flux distributions, Fi, oc i/"™". Pre- 
vious composite spectra from small numbers of AGN are shown in 
red (HST) and light blue (FUSE), with slopes as = 1.76 ± 0.12 
(Teller et al. 2002) and Os = O.bGtoH (Scott et al. 2004). Soft 
X-ray data are from ROSAT observations of AGN in the Lockman 
Hole (Hasinger 1994) normalized at i? = 1 keV (73.53 ryd) with 
spectral slope Fe oc £;-0-96±0.ii rpj^^ ionization potentials of H I 
and He II and ionization energies to produce key metal ions are 
shown as arrows. The metagalactic backgrounds are from Haardt 
& Madau (2001, 2011). 

HST/STIS to provide critically evaluated column densi- 
ties for iVjii and Nqyi, and their absorption-line frequen- 
cies, dAf/dz. The new O VI data were included, together 
with our fits to /ovi (Z/Zq), to yield more accurate val- 
ues of rtl^^^\ 

With the ten-fold increase in UV sensitivity through- 
put of Cosmic Origins Spectrograph on HST, we should 
be able to do even better. We hope to use HST/COS 
to obtain high-quality data (S/N > 30) to search for 
additional baryons in weak absorbers and constrain 
the predicted flattening in the column density distri- 
butions of H I and O VI. Scaled to column densities 
(10"'^'^ cm~^)iVi3, these lines have equivalent widths of 
(12.5 mA)Afi3 for O VI A1032 and (54.5 mA)iVi3 for Lya. 
These weak-absorber surveys will require HST/COS sen- 
sitivity to 4 mA equivalent widths, which is achievable 
at S/N = 30 toward many bright AGN background tar- 
gets. We can also use COS to obtain better detections 
and statistics for broad Lya absorbers (BLAs) and the 
Ne VIII doublet (AA770.4, 780.3). The Ne VIII Hues are 
potentially more reliable probes of hot, coUisionally ion- 
ized gas than O VI, since Ne VIII requires 207 eV to pro- 
duce and is likely to be less contaminated by photoioniza- 
tion. The lower solar neon abundance, (Ne/O)0 « 0.15, 
makes the Ne VIII lines weak, and redshifts z > 0.47 
are needed to shift them into the HST/COS band. The 
BLAs have considerable promise for WHIM probes, as 
they do not require corrections for metallicity. They do 
require determining the neutral fraction, /hi, through 
careful modeling of the gas temperature and ionization 
conditions. 

It would also be helpful to verify the claimed X-ray 
detections of O VII in the WHIM (Nicastro et al. 2005) 
which are not confirmed in other data (e.g., Kaastra et al. 




Fig. 9. — Compilation of current observational measurements of 
the low-redshift baryon census (Section 3.3). Slices of the pie-chart 
show baryons in collapsed form (galaxies, groups, clusters), in the 
circumgalactic medium (CGM) and intercluster medium (ICM) 
and in cold gas (H I and He I). Reservoirs include diffuse pho- 
toionized Lya forest and WHIM traced by O VI and broad Lya 
absorbers. Blended colors (BLAs and O Vl) have combined total of 
25 ± 8%, accounting for double-counting of WHIM at 10^ - 10^ K 
with detectable metal ions. The collapsed phases (galaxies, CGM, 
ICM, cold neutral gas) total 18 ± 4%. Formally, 29 ± 13% of 
the baryons remain unaccounted for. Our simulations (Figure 2) 
suggest that an additional 15% reside in X-ray absorbing gas at 
T > 10"^ K. Additional baryons may be found in weaker lines of 
low-column density O VI and Lya absorbers. Deeper spectroscopic 
UV and X-ray surveys are needed to resolve this issue. 

2006; Rasmussen et al. 2007) or in re-analysis of the 
same data (Yao et al. 2012). The most critical observa- 
tions for the WHIM census may require a next generation 
of X-ray spectrographs to measure the weak absorption 
fines of O VII A21.602, O VIII A18.969, and other He- 
like and H-like lines of abundant metals (C V, C VI, 
N VI, N VII). As discussed by Yao et al. (2012), this 
requires high-throughput spectrographs {E « 0.3 — 1.0 
keV) with energy resolution E/AE > 4000 sufficient to 
resolve O VII absorbers with mA equivalent width. For 
weak lines, the predicted O VII equivalent widths are 
Wx = (2.88 mA)(Aroyii/10^^ cm'^). 

Finally, the IGM simulations can be improved in sev- 
eral aspects, in order to use them as even more reli- 
able predictors of IGM parameters. In this paper, we 
computed individual (ccU-by-cell) values of metallicity Z 
and ionization fraction /ovi, together with their statis- 
tical variations and co-variance. Our simulations (Smith 
et al. 2011) were performed with somewhat larger box 
sizes (50ft.~^ Mpc) and on a larger 1024'^ grid, compared 
with previous studies of low-z IGM thermodynamics. 
Our O VI conclusions appear to be robust, as gauged 
by convergence tests and runs with different methods 
of injecting feedback from star formation into the grid. 
Post-processing these simulations allowed us to provide 
more accurate corrections for the product, {Z/ZQ)foYi, 
as a function of column density, Nqvi- These statistics 
can be performed for other key ions (O VII, Ne VIII). 
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We will extend our simulations through higher resolu- 
tion, addition of discrete ionizing sources, and radiative 
transfer. We also will explore the the correct mixture of 
coUisional ionization and photoionization in the WHIM, 
a project that requires understanding the implications of 
different feedback mechanisms for injecting mass, ther- 
mal energy, and metals into the CGM. How these metals 
mix and radiate likely determines the thermodynamics 
of the surrounding IGM. Because there still exist con- 
siderable differences between how simulations treat the 
thermal state of the IGM, we will continue to push to 
higher resolution to capture the low-mass galaxies which 
are important for metal production and mass injection. 
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NASA and AST07-07474 from NSF) at the University 
of Colorado Boulder. We thank Eric Hallman, Michele 
Trcnti, John Stockc and Evan Tilton for comments on 
the manuscript, and Mark Voit for discussions on hot 
gas in clusters and galaxy halos. 



APPENDIX A: Lya FOREST BARYON CONTENT AND 
ERROR ANALYSIS 

Here, we summarize the current status of the error 
budgets for the major contributors to uncertainty in the 
diffuse Lya forest baryon census. Many of the scalings 

enter the formula for O.'j^^^ as the square root, a weak de- 
pendence arising because recombination theory predicts 
that A^Hi oc njf for highly ionized absorbers. Thus, the 
total baryon mass (proportional to Uh) depends on a 

1/2 

weight factor N-^l (see equation [6]). From this formu- 
lation, we can conduct an error-propagation analysis, to 
arrive at the overall uncertainty on the Lya contribution 

to fib in equation (9). This quantity depends on five 
parameters, which we list below with assigned errors. 

1. Distribution of H I column densities. Dan- 
forth & Shull (2008) measured the distribution of 

H I column densities, dJ\f /dz cx iVj^j'^. Based on 650 
low-redshift Lya absorbers, they found a slope /? = 
1.73 ± 0.04 and a baryon content of 29 ± 4% of Qb, 
integrated over 12.5 < logA^ni < 16.5. Our new 
survey (Tilton et al. 2012) finds ^ = 1.68±0.03 and 
a contribution of 28 ± 4% from the Lya-forest plus 
partial Lyman-limit systems to the baryon fraction. 
We assume 15% uncertainty in the column-density 
distribution statistics, which are increasingly un- 
certain at logA^Hi < 13.0. Extending the N^^''' 
distribution from logN = 12.5 to 12.0 would add 
another 3% to the baryon fraction (from 28% to 
31%) if there is no change in slope. 

2. Hydrogen photoionization rate. This rate, 
Th, depends on the radiation field normalization 
at 1 ryd (/q) and spectral slope (as) between 1.0- 
1.5 ryd. Previous AGN composite spectra found 
indices ranging from = 1.76 ± 0.12 (Telfer et al. 



2001) to as = 0.56to;g (Scott et al. 2004). New 
HST/COS composite spectra of AGN find an index 
as = 1.59 ± 0.20 (Shull et al. 2012b), close to the 
Telfer et al. (2001) value for radio-quiet AGN. Our 
calculations also take into account the (1 + z)^ ** 
rise in Th{z) from z = to z = 0.7, a fit to the 
recent calculations of Haardt & Madau (2012). We 
estimate the joint error on Th as ±50% from mod- 
els of the metagalactic radiation field from quasars 
and galaxies (Shull et al. 1999; Haardt & Madau 
2012). 

3. Characteristic scale-length of absorbers. The 

characteristic impact parameter, pioo. at A^hi ~ 
10^^ cm~^ in units of 100 kpc, is inferred by di- 
rect and indirect means, including comparing "hits 
and misses" of Lya absorbers along nearby sight 
lines and the cumulative distributions of absorbers 
with nearest-neighbor galaxies (Stocke et al. 1995; 
Shull et al. 1998). The frequency of Lya absorp- 
tion lines per unit redshift also implies 200-300 
kpc absorber cross sections, when associated with 
the space density of galaxies down to luminosities 
L « 0.01 - 0.03L* (Shull et al. 1996; Stockc et al. 
2006; Prochaska et al. 2011). We adopt an uncer- 
tainty of ±50% on pioo- 

4. Electron temperature. Temperature (Tg) enters 
through the square root of the hydrogen recombi- 
nation rate coefficient, a^^ oc X'-0-726^ Models 
of IGM photoelectric heating (Donahue & Shull 
1991) predict a range from 5000 K to 30,000 K. 
For z < 0.4, with small expected variations from 
photoionization, we adopt an uncertainty of ±30%. 

5. Hubble constant. This parameter has been mea- 
sured as i?o = 72 ± 8 km s~^ Mpc~^ (Preedman 
et al. 2001) and Ho = 73.8 ± 2.4 km s^i Mpc^^ 
(Riess et al. 2011), using distance scales based on 
Cepheids in galaxies with Type la supernovae. To 
be conservative, we adopt an error of ±5%. 

From standard error-propagation formulae, we write 

the relative error on Qj^ as the quadrature sum of rela- 
tive errors on the five parameters, weighted by the square 
of the exponents (1, 0.5, or 0.363) as they appear in the 
scaling (see equation 9): 



\ N ) 



(0.5)^ 



P 



(0.363)" 



(Try 
T ) 



(10) 



This formula gives a relative error (crjj/O) = 0.40, so 

that we can express 0^"^V^^i*°*^ = 0.28 ±0.11. Most 
(77%) of the error budget comes from uncertainties in the 
ionizing radiation field (Th) and characteristic absorber 
size (pioo)- 
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APPENDIX B; CURRENT BARYON CENSUS 

The following paragraphs summarize our current knowl- 
edge of the baryon-content in various components and 
thermal phases of the IGM. 

• Photoionized Lya Absorbers. For this paper, 

/TTT\ 

we adopt fll = 28 ± 11% based on the results of 
Danforth & Shull (2008), our new survey (Tilton 
ct al. 2012), and the systematic uncertainties dis- 
cussed in Section 3.2. The mid-range distribution 
in column densities, 13.0 < log A^^hi < 14.5, is fairly 
well characterized, but the numbers of high-column 
absorbers are small. Their contribution to il^ re- 
mains uncertain owing to corrections for their size 
and neutral fraction. At the low end of the column- 
density distribution, there could be modest contri- 
butions to the baryon content from weaker Lya 
absorbers. In current surveys, their numbers are 
increasingly uncertain at logA^Hi < 13 (our inte- 
gration was down to 12.5). For a power-law dis- 
tribution with P = 1.7, extending the distribu- 
tion from logiV = 12.5 down to 12.0 would in- 
crease the baryon fraction by another 10% (from 
28% to 31%). For this paper (Figure 9), we adopt 

fi^ = 28 ± 11% based on the results of Danforth 
& Shun (2008), Tilton et al. (2012), and the sys- 
tematic uncertainties discussed in Section 3.2. 

• WHIM (O Vl-traced). Previous FUSE sur- 
veys of intergalactic O VI absorbers (Danforth & 
Shull 2005; Tripp et al. 2006) found lower Umits 

of 5% and 7%, respectively for the contribution 
of this gas to the baryon inventory. These sur- 
veys assumed metallicity Z = OAZq and ioniza- 
tion fraction /ovi = 0.2. In 2008, three O VI 
surveys with HST/STIS (Danforth & Shull 2008; 
Tripp et al. 2008; Thom & Chen 2008) probed 
to lower O VI column densities. Based on the 
survey of 83 O VI absorbers (Danforth & Shull 
2008) and our more recent analysis of 111 O VI 
absorbers (Tilton et al. 2012), we adopt a baryon 

fraction = 17 ± 4%. The factor-of-two in- 

crease arises primarily from our revised corrections, 
{Z/ZQ)foYi — 0.01, for metallicity and O VI ion- 
ization fraction (see Section 3.1). There may be 
some, as yet undetermined, overlap of photoionized 
O VI with the Lya forest. 

• WHIM (BLA-traced). Broad Lya absorbers 
(BLAs) were proposed (Richtcr et ai. 2004, 2006; 
Lehner et al. 2006, 2007) as repositories of a sub- 
stantial fraction of the low-redshift baryons. BLAs 
are defined as Lya absorbers with Doppler parame- 
ters & > 40 km s""'^, corresponding to temperatures 
T = {mHb^/2k) = (9.69 x W K)blo for pure ther- 
mal broadening with b = (40 km s~-'^)64o. Owing 
to the large abundance of hydrogen, a small neu- 
tral fraction (H I) remains detectable in Lya up 
to T 10^ K, although high-S/N is required to 
measure the broad, shallow absorption. As a re- 
sult, the surveys of BLAs differ considerably. In 
their survey of seven AGN sight lines, Lehner et al. 
(2007) found a BLA frequency of dAf/dz = 30 ±4 



for absorbers with 40 km s ^ < 6 < 150 km s ^ 
and logA^Hi > 13.2. They claimed that 20% of 
the baryons reside in BLAs. A more recent sur- 
vey (Danforth et al. 2010) came to different con- 
clusions. Surveying BLA candidates along seven 
AGN sight lines observed by HST/STIS, their 
BLA absorption-line frequency per unit redshift 
was dAf/dz = 18 ± 11, comparable to that of the 
O VI absorbers but 40% lower than that found by 
Lehnert et al. (2007). After correction for possible 
(20-40%) overlap between BLA and O VI (metal- 
bearing) absorbers, the corresponding baryon frac- 
tion is njiLA/^b = O-^'^-omt- For Figure 9, we 
adopt a value of 14 ± 7%, with an increased error 
reflecting the uncertain detection statistics. We 
apply their metallicity-based correction to obtain 
a blended total (approximately 1/4 overlap) of 
25 ± 8% for the O VI/BLA-traced WHIM over the 
temperature range 5 < logT < 6. 

• WHIM (X-ray absorber-traced). Our simula- 
tions (Figure 2) and those of Cen (2012) suggest 
that ^ 15% of the baryons could be contained in 
hotter WHIM, at T > 10^ K. Some of this gas 
could be detectable in X-ray absorption lines of 
trace metals (O VII, O VIII, Ne IX, NcX, etc.). 
However, most of these weak absorption lines are 
below the detection limits of current X-ray spec- 
trographs (Yao et al. 2012). The suggested absorp- 
tion systems are too few in number to provide good 
statistics, and many are unconfirmed. 

• Galaxies. Salucci & Persic (1999) found that 
galaxies contribute 7% of the baryons. More re- 
cent discussion by Fukugita & Peebles (2004) esti- 
mated 6%. In Figure 9, we assume that galaxies 
contribute 7 ± 2% of the baryons. 

• Groups and Clusters. The integrated clus- 
ter mass function of Bahcall & Cen (1993) was 
used by Fukugita et al. (1998) to find that the 
baryon contribution of clusters of galaxies consists 
^^^(stars) ^ 0.00155/1-1-5 and fll^'"''' = 0.003/1"!. 
Adjusting for h = 0.7, we find a cluster contri- 
bution of nl"^'' = 0.00308 or 6.8% of the baryons. 
Fukugita & Peebles (2004) revised the hot-baryon 

contribution to n^"^^ = 0.0018 ± 0.0007 or 4% of 

the baryons, owing to a redefinition of cluster mass 
(Reiprich & Bohringer 2002). For the total contri- 
bution of galaxy clusters, including their hot gas, 
we adopt a fraction 4.0 ± 1.5% of the baryons. 

• Cold H I (and He I) Gas. Zwaan et al. (2003) 
and Rosenberg & Schneider (2003) conducted blind 
H I surveys in the 21-cm line that probe the mass 
density in neutral atomic gas. The HIPASS survey 

(Zwaan ct al. 2003) foimd = (4.7±0.7) x lO"! 
Following the discussion of Fukugita & Peebles 
(2004), which augments the H I measurements by 
the expected accompanying cold He I and H2, we 
plot the total cold gas mass in Figure 9 as 1.7±0.4% 
of the baryons. 

• Circumgalactic Medium (CGM). X-ray 
spectra of AGN taken with both Chandra and 
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XMM/Newton detected strong O VII absorbers at 
z sa (see Bregman 2007; McKernan et al. 2005; 
Wang et al. 2005; Fang et al. 2006). The typical de- 
tected column densities, A^ovii ~ 10^^ cm~^, corre- 
spond to ionized hydrogen column densities ./Vhii ~ 
10^" cm~^ assuming a length-scale ~ 10 kpc and 
mean metallicity of 20% solar. Several arguments 
suggest that the Galactic O VII resides in a thick 
disk or low scale-length halo (5-10 kpc; see Yao 
& Wang 2005). Such a reservoir holds 10^ M© 
which is - 2% of the 5 x 10" M© in Milky Way 
baryons. Bregman (2007) suggested that a hot 
gaseous medium with mass ~ 10^" Mq might ex- 
tend throughout the Local Group. Such a reservoir 
is also ^2% of the 6 x 10^^ Mq baryonic mass of the 
Local Group, assuming total mass 5 x 10^^ Mq and 
12% baryon fraction (McGaugh et al. 2010). How- 
ever, the Milky Way halo and Local Group may 
not be typical of other star-forming galaxies. Re- 
cent COS studies of 42 galaxies (Tumlinson et al. 
2011) found that large (150 kpc) oxygen-rich ha- 
los of star-forming galaxies are major reservoirs of 
galactic metals. Additional CGM could extend to 
distances of 100-200 kpc from galaxies with higher 
specific star formation rates (Stocke et al. 2006). 
Savage et al. (2010) detected hot circumgalactic gas 
(logT = 5.8 — 6.2) in H I-free O VI absorption asso- 



ciated with a pair (perhaps small group) of galaxies 
at z w 0.167, at impact parameters p « 100 kpc. 
They note that the absence of H I absorption with 
broad O VI suggests a "rare but important class of 
low-z intergalactic medium absorbers" . 

Prochaska et al. (2011) took this idea further, sug- 
gesting that all the O VI arises in the "extended 
CGM of sub-L* galaxies". Their mass estimate 
associates the O VI absorption with H I around 
galaxies (0.1 < L/L* < L*) having constant hydro- 
gen column densities, Nh = 10^^ cm~^, over the 
full 200-300 kpc extent, as gauged by AGN- galaxy 
impact parameters. They estimated that these 
extended halos could contain a mass Mcgm ~ 
(3 X 10^° M0)(rcGM/3OO kpc)^, which would rep- 
resent ^ 50% of the baryon masses of present-day 
sub-L* galaxies. They further assert that gas at 
these large distances is gravitationally bound and 
virialized. We question the realism of constant- A^, 
virialized gas at such large distances from dwarf 
galaxies. Previous nearest-neighbor studies low-z 
O VI absorbers and galaxies (Stocke et al. 2006; 
Wakker & Savage 2009) found correlations with L* 
galaxies (at 800 kpc) and with dwarf (O.li*) galax- 
ies (at 200 kpc). For Figure 9, we adopt a CGM 
contribution of 5 ± 3%, recognizing that the CGM 
reservoir is still poorly understood. 
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